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Motivated by the recent finding of superconductivity in layered C0O2 compounds, we investigate 
' superconducting and magnetic instabilities of interacting electrons on the two-dimensional triangular 

, lattice. Using a one- loop renormalization group scheme for weak to moderate coupling strengths, 

f"^ ■ we find that for purely local interactions U > and small Fermi surfaces the renormalization 

£Nj ' group flow remains bounded down to very low scales and no superconducting or other instabilities 

> | , can be detected. Antiferromagnetic exchange interactions J generate a wide density region with a 

O i 1 d x 2_ y 2 +id xy -wa.ve superconducting instability similar to recent proposals for the strongly correlated 

I t-J model. For larger Fermi surface volumes the interactions flow to strong coupling also for purely 

local interactions U > 0. We find a singlet pairing instability in the vicinity of strong magnetic 
ordering tendencies at three wavevectors for the van Hove filling. 
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I. INTRODUCTION 



The triangular lattice has gained interest in modern condensed matter physics due to several reasons. First the 
triangular lattice continues to be important in the study of frustration effects in spin models, and is a promising 
candidate for the occurrence of spin-liquid phases 1 ' 2 . Then there are organic superconductors like the k-(BEDT-TTF) 
materials that may be described by an anisotropic triangular lattice and whose phase diagrams show some resemblance 
to those of the high-T c cuprates 3 . Very recently superconductivity was discovered at critical temperatures of 5K 4 
in Nao.35Co02 • I.3H2O. In this compound edge-sharing octahedra of C0O2 form two-dimensional layered triangular 
networks. A one-band Hubbard model on the two-dimensional triangular lattice may be the simplest approach to 
^vq ' describe the motion of the charge carriers in the layers. 

, In the cobalt oxides there are signs of electronic correlations like weak ferromagnetism and high temperature Curie 
behavior of the susceptibility 5 ' 6 . Therefore it is natural to ask whether the superconductivity in these materials 
, arises due to electron-electron interactions. The experimental findings have immediately spurred several theoretical 
suggestions 7-10 that RVB (resonating valence bond) superconductivity 11 might be realized in these materials. The 
RVB theories in the framework of the t- J-model consider the case of strong on-site repulsion on the Co sites. An 
alternative approach is to start from weak coupling. A very flexible way to investigate superconductivity arising from 
weak to moderate electron-electron interactions is to use renormalization group (RG) techniques. Modern functional 
RG studies have proved useful in the analysis of the Hubbard model on the two-dimensional square lattice 12-15 and 
give results that are in agreement with experiments and other theoretical approaches. In particular, on the square 
lattice for densities close to one electron per site, various strong coupling approaches 16,17 and weak coupling RG 
1 ■ calculations 12-15 agree on the symmetry of the superconducting order parameter. 

The Hubbard model on the triangular lattice has been analyzed using RG techniques by Tsai and Marston 18 . These 
authors were mainly concerned with the half-filled case and could not find any relevant superconducting tendencies 
for the isotropic triangular lattice. Vojta 19 and Dagotto studied superconducting instabilities due to paramagnon- 
exchange and found that the d xy channel dominates for most densities less than 1 electron/site. In addition they 
obtained an odd-frequency s-wave pairing as dominant channel close to half band filling. The fluctuation-exchange 
(FLEX) method was also applied to search for superconductivity on the anisotropic triangular lattice 20 . In this paper 
jJJ ■ we study two cases: 
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a) We analyze the parameter case that is, based on LDA (local density approximation) calculations for NaCo204 _ 
currently believed to be relevant for the C0O2 compound where superconductivity was found. NaCo204 has 1.5 
electrons per Co site in the A\ g -t2 g band close to the Fermi energy. For more details see the papers of Singh 21 and 
Wang et al. 9 . The superconducting samples (Nao.3sCo02 • I.3H2O) should have 0.15 electrons less per Co atom than 
in the NaCo204-layer, i.e. (n) w 1.35/Co site. If we want to describe the LDA band structure with a tight binding 
dispersion, the hopping amplitude t should be chosen negative. Then the Fermi surface (FS) is hole-like and closed 
around the origin, in agreement with photoemission data on NaCo2 4 22 . In our model calculation we choose t > 0. 
The case t < is obtained by interchanging particles and holes. Then the superconducting samples correspond to 
(n) ~ 0.65/site. We also study the half-filled band, (n) = 1/site, as this is interesting to make contact with other 
work. Furthermore we discuss the effect of exchange interactions between nearest neighbors. 
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b) The triangular lattice for t > at 3/4 band filling ((n) = 1.5/site) has a van Hove singularity in the density 
of states at the Fermi level, and the hexagon-shaped Fermi surface is nested. Thus we expect some kind of infrared 
instability of the Landau-Fermi liquid state, fc-space resolved A-patch RG may be the appropriate method for an 
unbiased description of the flow to strong coupling. 

The paper is organized as follows. In Sec. II we describe the ./V-patch renormalization group scheme and its 
implementation, in Sec. Ill we discuss case a) with small to medium band filling, and in Sec. IV we treat b), i.e. the 
density region near the van Hove filling. We conclude with some remarks on a possible relation to the cobalt oxides 
in Sec. V. 



II. THE MODEL AND THE METHOD 



The Hamiltonian of the Hubbard model on the 2D triangular lattice is 



H 



s Cj s + Cj s Ci >s 



Here t is the amplitude for hopping between nearest neighbors (ij) on the triangular lattice and U > is the on-site 
repulsion. We choose the reciprocal lattice to be spanned by the vectors (2tt, ±47r/V3) where we have set the lattice 
spacing to unity. Then the first Brillouin zone (BZ) is a regular hexagon (see Fig. 1). The non-interacting kinetic 
energy reads 



e(k) 



-2t 



cos k x + 2 cos — cos k y 

2 2 y 



where fi is the chemical potential regulating the particle density. The total bandwidth is 9t. For the cobalt oxides 
t < according to LDA calculations 21 . In this paper we will take t > 0. In order to treat the case t < we have to 
interchange particles and holes. 
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FIG. 1. Left: Brillouin zone and lines of constant band energy for the triangular lattice. The volumes encircled by these 
lines correspond to density n = 0.65, n = 1, n = 1.5 and n = 1.8 particles per site. Right: Density of states per spin and site. 
The chemical potential is chosen corresponding to n — 1.5 and the Fermi energy is at the van Hove density. 



Next we briefly describe the RG formalism that is used to obtain the effective interactions at low temperatures. We 
apply the so-called temperature-flow RG scheme introduced recently 15 in a iV-patch implementation that covers the 
full Fermi surface. Similar to the approaches in Refs. 12-14 the T-flow scheme is derived from an exact RG equation. 
However, no low energy cutoff is introduced. Instead of integrating out modes we use the temperature T itself as flow 
parameter and follow the evolution of the fermionic vertex functions as T is varied. As was argued recently 15 the 
procedure allows an unbiased comparison between AF and FM tendencies. In contrast with that, RG schemes with a 
flowing IR cutoff artificially suppress particle-hole excitations with small wave vectors, e.g. long wavelength density 
fluctuations. For a derivation and discussion of the T-flow RG scheme the reader is referred to Ref. 15. 

The RG formalism 23 produces a hierarchy of differential equations for the one-particle irreducible n-point vertex 

(n) 

functions, Tj, , as functions of the temperature. Integration of this system of equations gives the temperature-flow. 
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As initial condition we assume that at a high temperature To, the single-particle Green's function of the system is 
simply Go(ioj,k) = [iu — e(fc)] -1 and the interaction vertex is given by a local repulsion, U, or some other short 
range interaction. The initial conditions are justified if T is sufficiently large. We truncate the infinite system of 
equations by dropping all vertex functions T^?' with n > 4. In the present treatment we also neglect selfenergy 
corrections and the frequency dependence of the vertex functions. This restricts the scheme to one-loop equations for 
the spin-rotation invariant four-point vertex TjP- Starting with weak to moderate interactions, we follow the T-fiow 
of as T decreases. 

The four-point vertex is determined by a coupling function Vr(fci,fc2, k 3 ) (see Refs. 14, 15, 23). The initial 
condition for pure on-site repulsion is Vr (k\, fe, £3) = U . The numerical implementation follows the work of Zanchi 
and Schulz 12 and was already explained in Refs. 14 and 15. We define elongated phase space patches around straight 
lines from the origin to the boundaries of the BZ. Next we approximate Vr(fci, &2, fe) by a constant for all wave 
vectors in the same patch. We calculate the RG flow for the discrete subset of interaction vertices with each FS 
patch represented by a single wave vector. Most calculations were performed using 48 patches. We have checked 
the consistency of the discretization by comparing with up to 96 patches. Together with the flow of the interactions 
we calculate the flow of several static susceptibilities, as described below. This allows us to analyze which classes of 
coupling functions and which susceptibilities become important at low T. In many cases we observe a flow to strong 
coupling, i.e. at sufficiently low temperature some components of the coupling function Vr(fci, fo, Aft) become larger 
than the bandwidth. The approximations mentioned above fail when the couplings become too large. Therefore we 
stop the flow when the largest coupling exceeds a high value larger than the bandwidth, e.g. Vr.max = 18i. This 
defines a characteristic temperature T* of the flow to strong coupling. For superconducting and magnetic instabilities 
T* can be interpreted as estimate of the transition temperature if ordering is made possible, e.g. by including coupling 
in the third lattice direction. 

As in previous papers 15 we discuss the temperature flow at fixed chemical potential and not at fixed particle density. 
Since the relevant part of the flow only occurs at low temperatures we do not expect significant qualitative differences 
between the two possibilities. 



III. RESULTS AT HALF FILLING AND LESS 



First let us consider band fillings where the Fermi surface is well inside the first BZ with particle densities around half 
filling and less. We start the temperature-flow RG at higher initial temperatures ~ At with local on-site interactions 
U = 3.5t or At and follow the flow down to low scales. In Figs. 2 and 3 we display the behavior of the most repulsive 
and most attractive coupling constants. Although there is some flow, the coupling functions remain within the oder 
of the bandwidth 9t, even if we go down to very low scales. Above T <~ 1Q~H no tendency towards a flow to strong 
coupling, i.e. a divergence of the coupling constants, is observed. 

Quite generally, for a non-nested Fermi surface away from van Hove singularities, superconductivity is the only 
possible instability at low scales in a weakly coupled system. Kohn and Luttinger 24 gave an argument why supercon- 
ductivity arises even in the presence of repulsive interactions. Although in our case we do not find a singularity over 
a wide scale range, it is still interesting to investigate which pair-scattering channel is most attractive. This channel 
can in principle become unstable at even lower temperatures. However in real systems these weak instabilities with 
non-s-wave pairing symmetries may be wiped out easily by impurity scattering. 

For the two different densities in Figs. 2 and 3 we plot the pair scattering V(k, —k — > k' , —k') at low temperatures 
together with basis functions of the most attractive pairing channels. We have calculated the temperature-flow of 
the pairing susceptibilities with angular dependence sinm# and cos m6 for m = 1, 2, 3, . . . , N around the FS. For the 
half-filled case in Fig. 2 there are two harmonics that grow together more strongly than all other pairing channels. 
The dominant channels are d xy = sin 29 and d x 2_ y 2 = cos 29. These two <i-wave functions form the basis of Tg", 
which is 25 one of the two two-dimensional even-parity representations of the symmetry of the hexagonal lattice 26 . 
We expect that a possible superconducting pair amplitude will be proportional to a time-reversal symmetry-breaking 
superposition d x 2_ y 2 + id xy of these two basis functions as this maximizes the condensation energy. Lee and Feng 27 
applied the Gutzwiller approximation scheme and variational Monte Carlo to the Heisenbcrg model on the triangular 
lattice. Among the paired paramagnetic states they found the d x 2_ y 2 + id xy -w&ve state to have the lowest energy 28 . 
Remarkably this is the pairing symmetry as the one suggested - albeit with extremely low T c - by the weak coupling 
RG at small to moderate U/t. The same type of superconducting pairing was also found near half band filling by 
various RVB theories 7-10 . 

We note that in this perturbativc treatment of the on-site repulsion no precursors of a Mott transition at stronger 
U are visible. This is different from the Hubbard model on the square lattice close to half filling, where for a certain 
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parameter range the flow to strong coupling reveals tendencies towards incompressibility 




FIG. 2. Left: Flow of the most repulsive and most attractive coupling constant for /i — 0.85t ((n) = 1/site). The inset shows 
the FS in the first BZ. Right: The solid line shows the pair scattering V(k, — k — ► k', —k') at low T with the patch index k 
belonging to wavevector k around the 72 patches on the FS. The outgoing wavevector k' is fixed at point 3 (marked by the 
dashed line in the inset in the left plot). The dashed lines show the model pair scatterings V(9k,dk'=s) = — sin 29k sin 29^ and 
V{9k,9 k ' = 'i) = — cos 29k cos 26 '3 corresponding to the leading pairing susceptibilities with symmetries sin29k and cos9k- 
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FIG. 3. Left: Flow of the most repulsive and most attractive coupling constant for = — 0.7i ((n) « 0.65/site). The inset 
shows the FS in the first BZ. Right: The solid line shows the pair scattering V(k, —k — > k' , —k') at low T with the patch index 
k belonging to wavevector k around the 72 patches on the FS. The outgoing wavevector k' is fixed at point 1 (marked in the 
inset in the left plot). The dashed lines show the model pair scatterings V(0k,Ok'-i) = — sin 3#fc sin 3#i corresponding to the 
leading pairing susceptibilities with symmetries sin 36k ■ 

Next let us consider the case (n) sw 0.65/site that might, as explained in Sec. I, correspond to the superconducting 
cobalt oxide compounds. Here again no flow to strong coupling occurs down to very low temperatures T ~ 10 -5 i 
(see Fig. 3). The dominant pair scattering is in an odd-parity triplet pairing channel with angle dependence sin 38, 
corresponding to the Tl representation of Ref. 25. Comparing the basis function with the actual pair scattering 
rendered by the RG scheme, we notice that only two of the three repulsive portions of the sin 38 around the FS 
are reflected in the pair scattering. Thus even though this representation is the most attractive pairing channel, is 
does not fully match the effective interactions which mix in other symmetries at this scale. The second strongest 
components are d x 2_ y 2 and d xy , the first component is clearly visible for k' close to a maximum of d x 2_ y 2 in the right 
plot in Fig. 3. Again, a possible T c in this case may be extremely low. 

We repeat that for our parameters and purely local interactions < U < 4t we do not observe flows to strong 
coupling in the density range around and less than half filling. From the weak coupling perspective the reason simply 
is that the FS is relatively round and in the lack of nesting or van Hove singularities near the FS the one-loop 
corrections to the pair scattering do not generate any sufficiently attractive component in the pair scattering. Similar 
low energy scales for Kohn-Luttinger superconducting instabilities are found in the weakly repulsive Hubbard model 
on two-dimensional square lattice for a poorly nested Fermi surface away from the van Hove filling. 
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We note that with our RG scheme which neglects the frequency dependence of the interaction vertices, it is difficult 
to investigate odd-frequency pairing. This unconventional type of pair formation 29 was suggested by calculations by 
Vojta and Dagotto 19 for the Hubbard model on a triangular lattice for larger U/t. Although we cannot rule out this 
possibility for the time being, it is difficult to see how this type of instability should arise for the given rather generic 
situation unless drastic strong coupling effects are at work. 

As has become clear, on-site repulsion alone does not give rise to superconductivity at reasonable energy scales for 
these band fillings. Therefore we extend our analysis and consider exchange interactions between nearest neighbors, 

Hj = J Si ■ Sj 

(ij) 

with Si = \ Y^s s' ^s,s /C l s Ci,s< ■ In principle these interactions will be generated by virtual hopping processes for local 
repulsion U ^> t. Since we treat all interactions perturbatively, this regime is not accessible with our RG method. We 
will therefore just assume that the initial action contains such an exchange term and analyze the consequences. 

First let us consider the case U = and J = t at half band filling. Now the interactions flow to strong coupling 
at T w 0.02£, and an analysis of the coupling functions and susceptibilities shows that the divergence is a very clean 
Cooper instability with again dominant d x i_ y i + id xy character. This is shown in Fig. 4. In the upper right plot we 

can see that only the coupling functions Vt (ki ,^2,^3) with total incoming wavevector k\ + ki = grow rapidly towards 
low T. All other couplings remain smaller than the bandwidth. The pair scattering V"(fc,— fc — > k',—k') follows a 
— (cos 29k cos 28k' + sin 29k sin 29k' ) form when the angles 9k, 8k' are varied around the Fermi surface. This means 
that d x 2_ y 2 and id xy components, which belong to the same representation on the triangular lattice, get amplified in 
equivalent ways in a BCS gap equation with this pair scattering. 
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FIG. 4. Upper left plot: Flow of the most repulsive and most attractive coupling constant for /1 = — 0.85t ((n) ~ 1/site), 
U — and J = t. The inset shows the FS in the first BZ. Upper right: Snapshot of the coupling function Vr(fci , &2, £13) a ^ 
low T. The first outgoing wavevector £3 is fixed at point 4 and the wavevectors k\ and U2 vary around the Fermi surface. 
The diagonal features are the Cooper scattering processes with total incoming wavevector ki + ki = 0. The scale-bar indicates 
the strength of Vr(ki, ki, k$). Lower left: The symbols show the pair scattering V(k, —k — » k' , —k') at low T with the patch 
index k of wavevector k around the 48 patches on the FS. The outgoing wavevector k' is fixed at point 10 (circles in upper 
panel, see inset of upper left plot) and point 4 (squares in lower panel). The solid line shows the model pair scatterings 
V(9k,9y) oc — (sin 29k sin 26y + cos 29k cos 29k') corresponding to the leading pairing susceptibilities with symmetries sin2#fc 
and cos2#fc. Lower right: Temperature scale for the d x 2_ y 2 + id xy instability versus density for (7 = and J = t > 0. 
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The same picture occurs for general band fillings, ft turns out that at least for J = t, U = 0, the d x 2_ y 2 + id xy -wetve 
instability occurs over a wide density range from (n) = 0.6/site up to (n) = 1.8/site (see lower right plot of Fig. 4). 
The temperature scale for this instability is maximal near the van Hove filling. For J = t it gets smaller than 10 -5 £ for 
densities less than (n) = 0.6/site, i.e. close to the density range (n) w 0.65 relevant for the superconductivity in the 
cobalt oxides. Note however that this temperature scale depends exponentially on J and the density of states close 
to the Fermi level. We expect that an additional U > docs not change the leading instability for these densities. 
We have checked this for U = 2t, J = 0.5t. Of course, our weak coupling analysis is not capable of describing a Mott 
state at half filling for large U/t. 

It is interesting to observe the difference between the square lattice and the triangular lattice with nearest neighbor 
hoppings near half-filling. On the square lattice the flows to strong coupling and the suggested phases, mainly 
antifcrromagnetic order and d x - y 2-wave superconductivity come out of the weak coupling RG very similarly for pure 
on-site U or nearest neighbor exchange J. On the triangular lattice the Fermi surface is not nested and no particle-hole 
symmetry exists. The flows for pure U or pure J interactions differ strongly in the energy scales and, further away 
from half filling, also in the dominant pair scattering. 

Ferromagnetic exchange J = —t induces a weak tendency towards p-wave superconductivity of k+ik y -wave symme- 
try at 1/2 and 1/3 band filling. Around the van Hove filling we find a ferromagnetic instability. 

IV. AROUND THE VAN HOVE FILLING 

For chemical potential /j, = 2t or 3/4 band filling, the non- interacting Fermi surface forms a hexagon. Its corners 
coincide with the van Hove singularities at the center points of the 6 flat sides of the BZ. Therefore we may expect 
some kind of infrared instability giving rise to a flow to strong coupling of a certain class of coupling constants. In 
addition to the diverging density of states, the flat sides of the Fermi surface are perfectly nested, and this can cause 
a RG flow to strong coupling as well. The situation is reminiscent of the simple nearest-neighbor-hopping Fermi 
surface of the 2D square lattice at half filling, where two pairs of nested FS sides meet at the van Hove points at 
(±7r, ±7r). However on the triangular lattice we expect that the nesting between the three pairs of opposite flat sides 
is somewhat less effective. The three nesting wavevectors (after addition of a reciprocal lattice vector) point to the 
three inequivalent centers of the BZ sides. Thus they do not correspond to the same wavevector as in the square 
lattice case. Therefore the low energy phase space weight gets distributed on three different spots in the BZ. 

In Fig. 5 we summarize the flow at the van Hove filling for N — 96 patches. In the middle plot one observes several 
nesting features, i.e. bright, mainly vertical lines that correspond to fixed wave vectors transfer k 3 — fc 2 between first 
outgoing and second incoming particles (the spin indices of first incoming and first outgoing have to be the same). 
The most divergent processes belong to umklapp scattering processes between the vicinity of (but not directly at) 
the van Hove singularities at the points where the FS touches the BZ boundary. The ki and £3 corresponding to 
the most divergent processes are connected by dashed lines in the left plot of Fig. 5. These processes drive the spin 
susceptibility, and correspondingly the spin susceptibility seems to diverge at these wavevectors. In the right plot of 
Fig. 5 we show the q dependence of the static spin susceptibility along the line q x = 0, q y = . . . 2tt/\/3. Compared 
to other regions in g-space, it is enhanced along the whole line q x = 0, q y > when the temperature is decreased, 
but shoots up most rapidly very close to q y = 2n/V3 w 3.63 at the edge of the BZ. This point and the other two 
symmetry-related points where the spin susceptibility diverges with equal strength are marked by crosses in the left 
plot in Fig. 5. Since these wavevectors are equal to the three nesting vectors between the flat FS sides, we will 
denote them Qn- Due to the limited (f-space resolution of the iV-patch discretization we are not able to decide clearly 
whether the peak occurs directly at or somewhat away from the BZ boundary. Going from TV = 48 to N — 96 patches 
moves the peak a bit closer to the BZ edge, but for a clear answer not only the angular but also the radial dependence 
of the coupling constants must be considered. The susceptibility for ferromagnctism at q = grows as well towards 
T — > 0, but the combination of large density of states and nesting makes the instability at Qn far stronger. 

A divergence of the spin susceptibility at a given wave vector signals strong ordering tendencies at the given wave 
vector, either in the ground state or, if possible, at finite temperatures. In our case each of the Qm would correspond 
to antiferromagnetic ordering along two directions of the triangular lattice and ferromagnetic ordering along the third 
lattice direction. From our calculation we cannot tell whether the system will choose one of these possibilities, a 
certain superposition the three or whether ordering occurs at all. An analysis of the fluctuations around each of these 
ordcrings may give an answer to this question. 
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FIG. 5. Left: The first Brillouin zone of the triangular lattice and the Fermi surface (small circles for a N = 96-patch 
calculation) at the van Hove filling \i — 2t, (n) = 1.5/site. The broken lines indicate the wavevector transfers that grow most 
strongly towards low temperatures. The crosses denote the 'hot' wavevectors Qn near which the spin susceptibility grows most 
rapidly. Middle: Snapshot of the coupling function Vr(fci, £2, fe) a t l° w T. The first outgoing wavevector kz is fixed near a 
van Hove point at point 3 (diamond in left plot), and the wavevectors fci and hi vary around the Fermi surface. The most 
repulsive scattering occurs for hi = 78 (square in left plot). Right: Spin susceptibility Xb(q) a t low T along the line q x = 0, 
q y = 0. . . 27r/3. The enhancement is strongest near the nesting wavevector Qn with q y = 27r/v3 ~ 3.63 at the edge of the BZ. 



Next we decrease the particle density away from the van Hove filling. This reduces the density of states at the Fermi 
surface and cuts off the nesting. If the filling is sufficiently away from the van Hove filling we observe a saturation of 
the spin susceptibility at all wave vectors before the interactions get larger than twice the bandwidth, which is our 
criterion for stopping the flow and beyond which the approximations may fail. For U = At this happens roughly for 
densities ±0.03/site and further away from the van Hove filling of 1.5 electrons/site. Obviously, like in the RG studies 
of the Hubbard model on the 2D square lattice 12,13,15 , the precise boundaries of the spin density wave instability 
region depend somewhat on the scale where the susceptibilities are compared and on how far we trust the one-loop 
flow. 

The pairing channels are not cut off, and if the flow in the particle-hole channel above the saturation is still strong 
enough, it generates a sizable attractive component in the pair scattering. This component will flow in the particle- 
particle and eventually diverge in a Cooper instability at low temperatures. In Fig. 6 we plot the flow of the extremal 
coupling constants. We clearly observe the run-away flow of the coupling constants at T » 2.5 • 10 _5 t that is solely due 
to the Cooper pair scattering processes. All other couplings, including those in the spin density wave channel, saturate 
at higher temperatures, as can be seen from the long plateau in the flow of the most repulsive coupling constant before 
the Cooper processes take off. In the right plot of Fig. 6 we display the pair scattering around the Fermi surface. 
It has even parity in fc-space and roughly follows a sin 68 behavior, corresponding to the representation of D^h 
according to Ref. 25. The 12 nodes of the gap function are at the van Hove points and in the centers of the flat FS 
sides. The second strongest component in the pair scattering varies like sin 126* around the Fermi surface. The scale 
of the instability towards this unconventional i-wave state depends strongly on the vicinity to the van Hove filling. 
This is shown in Fig. 7. For larger on-site repulsion U, the density region with relevant values for T* becomes wider. 

As we have seen at the end of the last section, an antiferromagnetic exchange interaction J leads to a superconducting 
instability towards a time-reversal symmetry breaking state even at the van Hove density. Obviously there will be a 
change in the flow to strong coupling near the van Hove density when we change the initial interaction from U > 0, 
J = to U = 0, J > 0. For U = 2t and J = 0.5t, the instability is of d x 2_ y 2 + id xy -wa.ve type. Here we do not 
attempt to determine the boundary between the spin density wave instability in the first case and the superconducting 
instability in the latter. Note however that this may be an interesting question for an approach that is capable of 
describing how the exchange term is generated by the Hubbard interaction. 
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FIG. 6. Left: T-flow of most repulsive and most attractive coupling constant slightly away from the van Hove filling at 
fi = 1.94t and initial U = 3.5t. The inset shows the 96 points on the FS used in the TV-patch calculation. Right: pair scattering 
Vr(k, — k — > k' , —k') at low T with k' fixed at point 2 (marked with a diamond in the inset in the left plot) near but not at the 
van Hove points and k varying through the patches around the FS. 




FIG. 7. Left: Temperature scale T* for the flow to strong coupling versus chemical potential [i for initial U = At (solid line 
with circles) or U = 3.5t (dashed line). Around the van Hove filling (n = 2i) between the vertical broken lines (between densities 
1.48 and 1.53 per site) the leading instability for U = At is in the spin channel near wavevectors Qn, outside these lines the 
dominant instability is singlet superconductivity. Left: Flow of the leading spin (near Qn, dashed lines) and superconducting 
susceptibilities (solid lines) for /i = 1.98t (thick lines) and fj, — 1.96i (thin lines) for U = At. 



V. CONCLUSIONS 



We have analyzed the effective scattering processes at low temperatures of interacting electrons on the two- 
dimensional triangular lattice. We have used a A^-patch temperature-flow renormalization group scheme in the 
one-loop approximation that has been found useful 12,13,15 in the analysis of the weakly coupled Hubbard model on 
the 2D square lattice. Starting with a weak to moderate on-site repulsion U ~ 3 . . . At at high temperatures, we have 
investigated possible magnetic and superconducting instabilities when the temperature is lowered. For most fillings 
and purely local initial interactions, especially for densities less than wl.4 electrons per lattice site (for hopping pa- 
rameter t > 0), the interactions remain within the order of the bandwidth. Down to very low temperatures T ~ 10 _5 i 
no flows to strong coupling are observed. Replacing the pure on-site Hubbard interaction by an antiferromagnetic 
exchange interaction J, the RG flow leads to a superconducting instability towards a time-reversal symmetry breaking 
d x 2_ y 2 + id xy -w&ve state at sensible scales for a wide density range between (n) > 0.6/site and (n) > 1.8/site. This 
state has the same pairing symmetry as the superconducting RVB states found in recent theories for strong on-site 
repulsion 7-10 . It may also be obtained by ordinary mean-field theories that do not invoke spin-charge separation. 

Near the van Hove density (n) = 1.5/per site, we observe a flow to strong coupling also for purely local Hubbard 
interactions. The rapid growth of the interactions at low temperatures is driven by umklapp processes between the 
van Hove regions and nesting between the flat parts of the (nearly) hexagon-shaped Fermi surface. In the immediate 
vicinity of the van Hove filling the spin susceptibility near the wavevector q = (0, 27r/3) and two symmetry- related 
wavevectors diverges most strongly, signalling magnetic ordering tendencies at these wavevectors. It may be interesting 
to investigate whether the systems tends to order (if we allow coupling between layers) by choosing one wavevector 
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or by forming a superposition of the three. Slightly away from the van Hove filling (outside the range 1.5 ± 0.03/sitc 
for U = At), the flow in the spin density wave channel is cut off, and the dominant instability occurs in a singlet 
pairing channel with symmetry A(9) oc sin 60, having 12 nodes and 12 extrema when the angle moves around the 
Fermi surface. This superconducting instability occurs on both sides of the van Hove filling and the characteristic 
temperature T* falls of rapidly with increasing distance to the van Hove density. Since the superconducting state has 
nodes in the vicinity of the van Hove points, it may be susceptible to further instabilities at lower temperatures which 
try to remove the large density of states at zero energy. 

Turning to the reported superconductivity in the layered cobalt oxides, we have to keep in mind that it is far 
from settled that these materials can be described by a simple one-band Hubbard model on the triangular lattice. 
Neither is it clear that they can be treated using weak-coupling techniques. In fact estimates 21 of the bandwidth and 
Hubbard U for NaCo 2 04 seem to show that the material is strongly correlated. Nevertheless it may be helpful to 
know the weak coupling picture. Band structure calculations for NaCo204 by Singh 21 indicate that the main Fermi 
surface is hole-like and centered around the (0, 0) point. Angular resolved photoemission data support this picture 22 . 
Translated into the language of a tight-binding model this case corresponds to a negative hopping parameter t < for 
electrons. The superconducting samples may be at wl.35 electrons/Co site, or, with our sign convention t > that 
interchanges electrons and holes (n) = 0.65/site. This case is investigated in Section III. For a purely local Hubbard 
repulsion U > without additional exchange interactions no sizable superconducting tendencies can be detected with 
our weak coupling technique. Antiferromagnetic exchange interactions J > generate a superconducting instability 
towards a d x 2_ y 2 + id xy -wa,ve state over a wide density range. In this respect, weak coupling RG and RVB theories 
for the strongly correlated t-J model 7-10 give similar predictions at least as far as the symmetry of the pairing is 
concerned. Thus the time-reversal symmetry breaking d + id' state may be the primary candidate for experimental 
checks. This state has a fully gapped single-particle excitation spectrum. The phase structure of the order parameter 
may give rise to observable effects 30 such as subgap Andreev bound states at surfaces similar to those in p-wave 
superconductors 31 ' 32 . 

The main purpose of the present work was to investigate possible superconducting states on the triangular lattice at 
weak to moderate interactions. Further work is needed in order to clarify the relevance of the weak coupling physics 
to real materials like the cobalt oxides. 

I thank P.A. Lee, V. Liu, O. Motrunich, T.M. Rice, T. Senthil, and A. Vishwanath for helpful discussions. The 
Deutsche Forschungsgemeinschaft (DFG) is acknowledged for financial support. 
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